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Abstract: We investigate the effects of the accretion of phantom energy with 
non-zero bulk viscosity onto a Schwarzschild black hole and show that black holes 
accreting viscous phantom energy will lose mass rapidly compared to the non- viscous 
case. When matter is incorporated along with the phantom energy, the black holes 
meet with the same fate as bulk viscous forces dominate matter accretion. If the 
phantom energy has large bulk viscosity, then the mass of the black hole will reduce 
faster than in the small viscosity case. 
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I. INTRODUCTION 

Observations of WMAP [H [21 G] and supernova of type la data pf] have 
revealed that our Universe is filled with an exotic dark energy apart from 
dark matter. The nature and composition of this energy is still an open 
problem but its dynamics is well understood i.e it causes an approximately 
exponential expansion of the Universe (see [5] for recent reviews on dark en- 
ergy). Astrophysical data suggest that about two thirds of the critical energy 
density is stored in the dark energy component. For the equation of state 
(EoS) parameter u < — 1, the fluid is called phantom energy (PE). Observa- 
tions show that u is constrained in the range —1.38 < u < —0.82 thus 
providing evidence of phantom energy in the Universe. The PE violates all 
the energy conditions in all forms (weak, null, strong or dominant). The 
phantom energy can cause some peculiar phenomena e.g. the existence of 
wormholes (TJ, Ej, infinite expansion of the Universe in a finite time causing 
a Big Rip (BR) and the destruction of all gravitationally bound structures 
including black holes P, [TUJ, dH [T2]. In particular, black holes will contin- 
uously lose mass and disappear near the BR (see [EH [H] for the opposite 
viewpoint). 

Dark energy with bulk viscosity has a peculiar property to cause acceler- 
ated expansion of phantom type in the late evolution of the universe [To] . 
It can also alleviate several cosmological puzzles like cosmic age problem 
[16J, coincidence problem |T7j and phantom crossing [TBJ. We will consider 
phantom energy as an imperfect fluid, implying that the PE could contain 
non-zero bulk and shear viscosities [19] . The bulk viscosities are negligible 



for non-relativistic and ultra-relativistic fluids but are important for the in- 
termediate cases. In viscous cosmology, shear viscosities arise in relation 
to space anisotropy while the bulk viscosity accounts for the space isotropy 



T51 120]. Generally, shear viscosities are ignored (as the CMB does not indi- 
cate significant anisotropics) and only bulk viscosities are taken into account 
for the fluids in the cosmological context. Moreover, bulk viscosity related 
to a grand unified theory phase transition may lead to an explanation of the 
accelerated cosmic expansion |2Tj . 

Babichev et al [TO] studied the effects of the accretion of phantom energy 
onto a Schwarzschild black hole taking PE to be a perfect fluid. As a first 
approximation, the bulk viscosity can be ignored, but to get a better picture 
we need to incorporate it into the phantom fluid. We have adopted the 



procedure of [TU|, E2J for our calculations. 



The plan of the paper is as follows: in the next section we review viscous 
cosmology; in section three we discuss the relativistic model of accretion onto 
a black hole; in the subsequent section we use results from viscous cosmology 
for the accretion model; next we give two examples to illustrate the accretion 
process with a constant and power law viscosity. In section six we study black 
hole evolution in the presence of matter and viscous phantom energy. Finally 
we conclude the paper with a brief discussion of our results. 

II. BULK- VISCOUS COSMOLOGY 



We assume the background spacetime to be homogeneous, isotropic and 
spatially flat (k = 0) and described by the Friedmann-Robertson- Walker 



(FRW) metric given by 

ds 2 = -dt 2 + a 2 {t)[dr 2 + r 2 {d6 2 + sin 2 0# 2 )], (1) 

where a(t) is the scale factor. We also assume that the spacetime is filled 
with only one component fluid i.e. the viscous phantom energy of energy 
density p (however, in section six, we shall incorporate matter along with 
phantom energy). The Einstein field equations for the FRW-metric (in the 
units c = 1 = 8ttG) are 

"M-) 2 = ^, (2) 



and 

- = H + H 2 = -±(p + Sp), (3) 
a 6 

where H is the Hubble parameter, p is the effective pressure containing the 
isotropic pressure p pe and the bulk viscous pressure p v - ls , given by 

P = Ppe + Pvis- (4) 

Here p = p pc + p v [ s and p v [ s = — where is the velocity four vector 
and £ = £(p y \s,t) is the bulk viscosity of the fluid [23] . Eq. (4) shows that 
negative pressure due to viscosity contributes in the effective pressure which 
cause accelerated expansion. In the FRW model, the expression ufi = 3d/ a 
holds. Also, £ is generally taken to be positive in order to avoid the violation 
of second law of thermodynamics [24 . 



The energy conservation equation is 

p + 3H(p + p) = 0. (5) 
Assume that the viscous fluid equation of state (EoS) is 

p = up = ( 7 - (6) 



Note that if 7 = (or u = — 1), Eq. (6) represents the EoS for cosmological 
constant. Furthermore if 7 < 0, it represents phantom energy. In general, 
for normal matter 1 < 7 < 2. 

Using Eqs. (2) - (6), we get the equation governing the evolution of H(t) 
for a given £ as 

2H + 3jH 2 - 3£H = 0. (7) 
On integration, Eq. (7) gives 

exp{ljgt)dt} 
C+|7exp{|/£(^} : 

where C is a constant of integration. Note that Eq. (8) can further be solved 

to get the evolution of a(t) as 



m = „ .T''^;!,, (8) 



a(t) = d(c + j exp j Z(t)dt}dt)°\ (9) 



where D is a constant of integration. Thus for a given value of £ we can 
obtain expressions of a(t), pit) and p(t) from the system of Eqs. (5) - (9). 

III. ACCRETION ONTO BLACK HOLE 

In the background of FRW spacetime, we consider, as an approximation, 
a gravitationally isolated Schwarzschild black hole (BH) of mass M whose 
metric is specified by the line element: 

ds 2 = - (l - ^jdt 2 + (l - ^) ~'dr 2 + r 2 (d9 2 + sin 2 6dtp 2 ). (10) 

The background spacetime is assumed to contain one test fluid, namely the 
phantom energy with non- vanishing bulk viscous stress p y [ s . The fluid is 
assumed to fall onto the BH horizon in the radial direction only which is 



in conformity with the spherical symmetry of the BH. Thus, the velocity 
four vector of the phantom fluid is w M = (V(r), u r (r), 0, 0) which satisfies the 



normalization condition u^u, 



1. This phantom fluid is specified by the 



stress energy tensor for a viscous fluid P2 



= {p + p)u»u v + pg^. 



Using the energy momentum conservation for T^ v ', we get 



(11) 



ur 2 M' 2 {p + p)\jl 



M 
47rr 



(12) 



where u r = u = dr/ds is the radial component of the velocity four vector and 
C\ is a constant of integration. The second constant of motion is obtained 
by contracting the velocity four vector of the phantom fluid with the stress 
energy tensor u^T.^f = 0, which gives 

ph 



ur M~ exp 



dp' 



ft + p(p') 



= -A, 



(13) 



where A is a constant of integration. Also ph and p^ are the energy densities 
of the phantom fluid at the horizon of the BH, and at infinity respectively. 
From Eqs. (12) and (13) we have 



M 

47TT 



Ph 



+ u 2 exp 



dp' 



ft + p{p' 



Co 



(14) 



with C2 = —C\/A = poo + p(poo)- In order to calculate the rate of change 
of mass of black hole M, we integrate the flux of the bulk viscous phantom 
fluid over the entire BH horizon to get 



M = <b T[dS. 



(15) 



Here T[ determines the energy momentum flux in the radial direction only 
and dS = ^—gdOdip is the infinitesimal surface element of the BH horizon. 
Using Eqs. (12) - (15), we get 

dM AM 2 . 

which clearly demonstrates the vanishing mass of the black hole if p + p < 0. 
Integration of Eq. (16) leads to 

M = Mo(l- t -)~\ (17) 

where Mq is the initial mass of the black hole and modified characteristic 
accretion time scale r _1 = [^|^{(ppe + Ppe) ~~ ^ m (^))]' a o being the value 
of the scale factor at time to. Note that during the integration of (16), we 
assumed p pc and p pc to be constants. In the coming subsections, we shall 
take these as time dependent entities. 

IV. ACCRETION OF VISCOUS PHANTOM FLUID 

We now study the BH mass evolution in two special cases: (a) constant 
viscosity; and (b) power law viscosity. 

A. Constant bulk viscosity 

For constant viscosity £ = £ , the evolution of a(t) is determined by using 
Eq. (9). It gives 

a(t) = a H , (18) 

where 

'3*6 



Bit) ee exp - 1. (19) 
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Using Eqs. (5), (6) and (8) the density evolution is given by 

p exp (3£ t) 



P(t) = 



1 + 



jHoBjt) 



(20) 



Here p Q = 3i7 2 . Further, for 7 < the BR singularity occurs in a finite time 
at 

Finally, the BH mass evolution is determined by solving Eq. (16) and (20) 
to get 



M = M 



where 



3tjo 



A = C + (-l + e— ) 1 H 



(22) 



(23) 



This mass is displayed for different values of viscosity at different times in 
Table 1. 





6 = 10- 17 


6 = 10- 18 


6 = 10- 19 


U = 10- 20 


ti = 10 7 


3.43427 x 10~ 4 


2.44662 x 10~ 3 


6.31285 x 10~ 3 


7.49184 x 10~ 3 


h = 10 10 


3.43544 x 10~ 7 


2.45261 x 10~ 6 


6.35248 x 10~ 6 


7.55357 x 10~ 6 


h = 10 13 


3.43516 x 10~ 10 


2.45258 x 10~ 9 


6.35247 x 10~ 9 


7.55358 x 10~ 9 


u = 10 17 


1.23994 x 10~ 14 


2.10182 x 10~ 13 


5.86096 x KT 13 


7.01997 x KT 13 



Table 1. The mass ratio M/Mq of black hole for different choices of 
constant viscosity £ . The initial mass is, throughout, taken to be 5OM or 
10 32 kg. 

It is apparent from Table. 1 that for a fixed viscosity, the mass ratio de- 
creases with time implying that mass of black hole is decreasing for an initial 
mass. Similarly, at any given time, the mass ratio also decreases with the 
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increase in viscosity. Thus the greater the value of viscosity parameter, the 
greater would be its effects on the BH mass. 



B. Power law viscosity 

If the viscosity has power law dependence upon density i.e. £ = apy is , 
where a and s are constant parameters, it has been shown [ 1251 127] that it 



yields cosmologies with a BR if y/Sa > 7 and s = 1/2. Thus we take £ = ap? 
as a special case. Then the scale factor evolves as 

a(t) = oo(l- J)^. (24) 
The density of phantom fluid evolves as 

^ s^-W 1 '^ (25) 

or in terms of critical density p CT as 

p(t) = pJl--Y 2 . (26) 



r 

The corresponding BR time r is given by 

2 

r 



H~\ (27) 



3(v3a: — 7) 

Finally, the mass evolution of BH is determined by using Eq. (16) and (25) 
is 

M = M \l+ A ^ ' - I" 1 . (28) 

L 4 7 r(v / 3a - 7) r(r - t) J V ' 

Note that when a = 0, this case reduces to that of Babichev et al [TO] . The 
mass in (28) in displayed for different values of EoS parameter 7 at different 
times in Table 2 and displayed graphically in Figure 1. As shown, the mass 



10 



decreases gradually with the decrease in the EoS parameter 7. Note that we 
have not graphically displayed the mass for different viscosities given in Table 
1 because the variation is not significantly different for most time scales. 





71 = -1 x 10- 1 


72 = -2 x nr 1 


73 = -3 x HT 1 


74 = -4x HT 1 


h = 10 10 


4.71915 x 1CT 5 


2.35963 x 10~ 5 


1.5731 x KT 5 


1.17983 x 10~ 5 


H = 10 13 


4.79136 x 10~ 8 


2.35968 x 1CT 8 


1.57312 x 10~ 8 


1.17984 x 10~ 8 


h = 10 17 


4.66492 x 1CT 12 


2.30523 x 10~ 12 


1.51867 x nr 12 


1.12539 x nr 12 


U = 10 20 


4.97349 x 10~ 14 


5.20946 x 10~ 14 


5.28811 x nr 14 


5.32744 x 10~ 14 



Table 2. The mass ratio M/Mq of black hole for different choices of 



equation of state. The initial mass is 5OM or 10 kg. 



V. EXAMPLES 

We now solve examples to demonstrate the accretion of viscous phantom 
energy onto a BH. The formalism is adapted from [TP] , 

A. Viscous linear EoS 

We choose the viscous linear EoS, p = cup pe — SH^ a with u < — 1. The 
ratio of the number densities of phantom fluid particles at the horizon and 
at infinity is given by 



n(p 



pe\ 

OO ; 



(29) 



.pg(l + w)-3&fT 

The critical points of accretion (the point where the speed of fluid flow be- 
comes equal to the speed of sound i.e. it 2 = c 2 ) are given by 



2 UJ 1 + 3u 



M * = TT3^ ; x ' = — ■ (30) 
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The constant A appearing in Eq. (16) is determined to be 

_ \i + M - 

A ~ 4|of/2 ' 

Notice that the constant A is the same as for the non- viscous case [10J. Also, 
the density of phantom energy at the horizon is given by 

p1 =t^ + \a) I^-tW- (32) 

Moreover, the speed of flow at the horizon is 

«* = -(£)'*'■ (33) 

The speed is negative as it is directed towards the BH. Also, the characteristic 
evolution time scale of the BH is given by 

Finally, substituting Eq. (34) in (17) we get the mass evolution of a BH in 
bulk viscous cosmology 

M = M„ [l - %M„t^±M 2 " {^e (1 + ul) _ 3|o ln Q jj - 1 (35) 

Since is unknown for our purpose, we have not evaluated M for different 
times numerically for tabular and graphical presentation. 

B. Viscous non-linear EoS 

We here choose the EoS, p = up pc — 3H^(p v i s ) with u < — 1, where £(p P e) = 
«/9p e with a and s are constants. The ratio of number densities is given by 



n(pm) \pooJ V + cj) - 3Hap 6 h 



(36) 
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The constant A appearing in Eq. (16) is determined to be 

A ( Ph\ (»-i)(i+*) / pooQ- + v) ~ SHap 3 ^ ^ 
\ Poo ) V p h (l+u) -3Hap s h 

The speed of flow at the horizon becomes 



(37) 



2s 



Ph \ 



V Ph (l + u) -3Hapi J 



, Poo / v p h (l + u) - 3Hap s h 

The critical points of accretion are given by 

2 u — 3sap s h l 1 + 3(cj — 2>sap s h ~ l ) 



(38) 



* 1 + 3(a; - 3sa/op) 2(a; - Ssapf 1 ) 

The characteristic evolution time scale r is given by 

t= 4ttMo — — — — r — — -) \ Poo 1+cj)-3— — In 

L Voo' V ^(1 + u;) - 3Hap s h J L t 

Finally, using Eqs. (40) in (17), the BH mass evolution is given by 



(39) 



a 



a 
(40) 



M = M 



1 - 47TM t 



ph \ (.-1)(1-U0 /Poo(l + CJ) - 3/^0!^^ 3 



V P/l (l+o;)-3#a^ 
x^( 1+ .)-3^.ne}; 



(41) 



As before, /Oqq is unknown, but further p^ is also unknown. As such, we again 
do not provide a tabular or graphical presentation. 



VI. BLACK HOLES ACCRETING BOTH MATTER AND VISCOUS 

PHANTOM FLUID 



We now consider a two component fluid, the viscous dark energy and mat- 
ter. The matter part may be composed of both baryonic and non-baryonic 
matter. It is taken to be a perfect fluid while the PE is taken as a bulk viscous 
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fluid. The effective pressure is represented by Eq. (4). The corresponding 
Einstein field equations (EFE) for the two component fluid become: 

1 



R, 



-a = T + T m 



The stress-energy tensor representing the two component fluid is given by 



(42) 



T"" = (p + p + p m )u'V + pg 



(43) 



Here p m is the energy density of the pressureless matter. Energy conservation 
holds independently for both fluids: 



p + 3H(p + p) = 0, 



Pm + SHp m = 0. 



(44) 
(45) 



Integrating Eq. (45), we have 



_3 

Pm Pmo® j 



(46) 



where p mQ = p m (to). Similarly, integrating Eq. (44) leads to 



p = p T 



3/2 



(47) 



where S is a constant and K is given by 

3>/3& 



X = 



yfP~ni Q 

Thus the total energy density of the two component fluid is given by [25] 

p = p + pm = pm Q a~ 3 + y fl3/2 ) • 
Using Eqs. (45) in (16) the evolution of black hole mass is given by 

-2 



(48) 



(49) 



M = M, 



1 - AnAM n 



9a 3 



(3S + 4^ /2 i^) 



^i^ln(^--f,(3- + 4a 3 /^) 
#(t)l 9 \a J 9a 3 v 

n -1 



(50) 
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where the scale factor a(t) evolves as 



a(t) 



-3 
.K 



( e f y/pm Q /it+D 1 _ ^ 



1 2 /3 



(51) 



and D\ is the constant of integration determined by choosing t = to get 



As pointed out in the next section, we cannot correctly discuss a BR scenario. 
However we can take a spacetime approximating it sufficiently earlier than 
the BR. We can than see its asymptotic behavior, when the scale factor 
shoots to infinity, the three terms in Eq. (50) will contribute significantly in 
the BH mass evolution. The mass will decrease by the accretion of PE (7 < 0) 
due to its strong negative pressure and is manifested in Eq. (50). Notice that 
the final expression for BH mass depends only on the initial matter density 
p mo in addition to constant bulk viscosity £ . The corresponding behavior of 
BH mass evolution is shown in Figures 2 and 3 for different values of model 
parameters. Thus for a shift of parameter 7 by 2, yields in the decline of mass 
ratio by a factor of 2. The decline in the mass of the BH is observed with 
time showing that phantom energy accretion will be dominant over matter 
accretion. 



We have analyzed the accretion of bulk viscous phantom energy onto a 
BH. The modeling is based on the relativistic model of accretion for compact 
objects. The viscosity effects in cosmology are used to give an alternative to 
cosmic accelerated expansion other then dark energy and quintessence. The 




(52) 



VII. CONCLUSION AND DISCUSSION 
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evolution of BHs in such a Universe accreting viscous phantom energy would 
result in a gradual decrease in mass. This gradual decline would be faster 
than the non-viscous case [TO] due to additional terms containing viscosities 
coupled with mass. Lastly, it is shown that BHs accreting both matter and 
viscous PE will also meet with the same fate as the viscous forces dominate 
over the matter component for sufficiently large scale factor a(t). 

From this analysis, we can draw the conclusion that PE containing viscous 
stresses can play a significant role in the BH mass evolution if the viscosity 
is sufficiently high for an appropriate EoS. Though the viscous stresses are 
negligibly small O(10~ 8 Nsm~ 2 ) at the local scale of space and time they can 
play a significant role in time scales of ~ Gyrs. The higher the viscosity of the 
phantom fluid, the sharper the decrease in the BH mass. BHs of all masses, 
ranging from the solar mass to the intermediate mass to the supermassive, 
will all meet the same fate. 

As an extension to this problem, it is interesting to study the accretion of 
the phantom fluid onto primordial BHs that had formed due to initial density 
fluctuations in the primordial plasma. The mini-primordial BHs evaporating 
now via Hawking radiation would have a different initial mass and hence 
abundance than the standard scenario expects. This work is reported in a 
separate paper [28] . 

Notice that we have used the Friedmann model which is represented by an 
asymptotically curved spacetime and at the same time the Schwarzschild 
black hole, which is asymptotically flat. This may seem contradictory. 
Schwarzschild black hole has been dealt with in the context of closed Fried- 
mann cosmology [212, EDI EI] • Any global problem in approximating the full 
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situation by a Schwarzschild black hole inserted into Friedmann model arise 
near the big bang or the big crunch, defined in terms of the york time [32] as 
shown elsewhere [33] , the effect will be at extremely late times in terms of the 
usual time parameter. More complete analysis of the asymptotic behavior 
near a singularity is also available [33] , as such if we stay near to a singularity 
in spacetime, the approximation will be extremely good. Consequently our 
analysis will be satisfactory for black holes formed well after the big bang 
greater then 10~ 40 s and of the Big Rip (presumably much more before 10~ 40 s 
the rip). It is clear that we are unable to say whether there would/would not 
be a Big Rip as our analysis excludes it. 
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FIG. 1: For an initial mass of black hole Mq = 10 32 kg, the evolution of the mass parameter 
m = M/Mq — 1 is plotted against the logarithmic time with a = 10~ 5 and tn = 10 17 s. 




2x 10 16 4x 10 16 6x 10 16 8x 10 16 lx ID 1 ' 



FIG. 2: For an initial mass of black hole Mo = 10 kg, the evolution of m is plotted against the time 
parameter t with A = 1/3, E = 3, £ G = lO'^fcym - ^ -1 and 7 = -1CT 1 while H w 2.33 x 10" 18 m. 
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FIG. 3: For an initial mass of black hole Mq = 10 32 kg, the evolution of m is plotted against 
the time parameter t with A = 1/3, S = 3, £ c = 10 _16 &;<7m _1 s _1 and 7 = —2 x 10" 1 while 
H w 2.33 x 10- 18 m. 



